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ABSTRACT 

Aims. Molecular hydrogen is the most abundant molecule in the Universe. It is thought that a large portion of H 2 forms by association 
of hydrogen atoms to polycyclic aromatic hydrocarbons (PAHs). We model the influence of PAHs on total H 2 formation rates in 
photodissociation regions (PDRs) and assess the effect of these formation rates on the total cloud structure. 

Methods. We set up a chemical kinetic model at steady state in a PDR environment and included radiative transfer to calculate 
the chemistry at different depths in the PDR. This model includes known dust grain chemistry for the formation of H 2 and a H 2 
formation mechanism on PAHs. Since H 2 formation on PAHs is impeded by thermal barriers, this pathway is only efficient at higher 
temperatures (T > 200 K). At these temperatures the conventional route of H 2 formation via H atoms physisorbed on dust grains is no 
longer feasible, so the PAH mechanism enlarges the region where H 2 formation is possible. 

Results. We find that PAHs have a significant influence on the structure of PDRs. The extinction at which the transition from atomic 
to molecular hydrogen occurs strongly depends on the presence of PAHs, especially for PDRs with a strong external radiation field. 
A sharp spatial transition between fully dehydrogenated PAHs on the outside of the cloud and normally hydrogenated PAHs on the 
inside is found. As a proof of concept, we use coronene to show that H 2 forms very efficiently on PAHs, and that this process can 
reproduce the high H 2 formation rates derived in several PDRs. 

Key words, astrochemistry - ISM: molecules - ISM: photon-dominated region 


1. Introduction 


Molec ular hydrogen is the most abundant molecule in the Uni¬ 
verse ( ^ull & Beckwithl 119821 : iBlack & van DishoeckI 119871 : 
Draine & Bertoldil ' 19961) and the main constituent of regions 
where stars are forming. Many ion-molecule and neutral-neutral 
reactions are driven by the presence of H 2 . Specifically, 
H 2 facilitates the formation of molecules such as CO, OH , 
H 2 O, HCO^, HCN, and Hj dOppenhei mer & D algar^ 1^[74|: 
van Dishoeck & BlackI ll986l:^Tielens & Hollenbachl 11985 ). 


These species, as H 2 itself, can be important coolants and 
provide usef ul diagnostics of (irradiated) interstellar clouds 
and shocks (Neufeld & KaufmanI 119931 : iNeufeld et al.l I1995L 

it is crucial to understand 
10 K) to high (~ 10^ K) 


Le Bourlot et al.l 19991) . As such. 


the formation of H 2 from low (' 


Since gas phase routes to form H 2 were found to be in¬ 
efficient, dust gr ains are recognized as the favoured habitat 
for H 2 formation (lOort & van de Hu'isill946l : iGould & Salneteil 
Il963l) . These dust grains can adsorb hydrogen atoms on 
their surface with either a Van d e r Waals bond (phys isorp- 
tion, iPirronello et al.l 1 19971 119991: iPerrv & Pricel 120031) or a 
covalent chemical bond (chemisorpt ion, IZecho et al.l l2002t 
iHornekaer et al.l[20061: lMennell^l20()6l) . The typical physisorp- 
tion energy of a hydrogen atom is 700 - 800 K (60 - 70 
meV), depending on the dust composition, whereas the typi- 
cal chemisorption energy is 7000 - 25000 K (0.6 0 - 2.15 eV) 
(iDulev & William^ I I 993 I ICazaux & Spaansll2009t) . Because of 
their lower binding energies, physisorbed H atoms can easily mi¬ 


grate to another site (physisorbed or chemisorbed), or can evap¬ 
orate back into the gas phase. If a migrating atom encounters 
another adsorbed H atom, these two atoms can react to form 
H 2 , the Langmuir-Hinshelwood mechanism (LH). H 2 can also 
be formed when a gas phase atom adsorbs onto an already occu¬ 
pied site (Eley-Rideal mechanism, ER). 

For both the LH and the ER mechanism to work, hydro¬ 
gen atoms must be able to adsorb and stick to the surface of 
the dust grain. As the grain temperature rises, physisorbed H 
atoms start to evaporate and only chemisorbed H atoms are left 
on the grain surface. This decreases the H 2 formation rate, since 
the LH mechanism becomes inefficient and H 2 formation re¬ 
lies on the ER mechanism involving chemisorbed H atoms. The 
temperature at which the change of efficiency occurs depends 
on the physisorption energy, which depends on the composi¬ 
tion of the dust grain. Typically, physisorbed H atoms desorb 
at grain temperatures of 10 - 20 K dHollenbach & S alpeteJ 19711 : 


ICuppen & Herbslj[2005l:lc azaux & Spaansll2009l). 


Howeyer, high H 2 formation rates haye been obseryed 
in PDRs with high dust and gas temperatures, such that 
the obseryed rates cannot be explai n ed by dust grains alone 
(iHabart et al.ll2004lAllers et al.l[200l . lAUe rs et al. 1 (1200-51) mod¬ 
elled the structure of the Orion Bar PDR and matched it with 
obseryations of different H 2 royibrational lines. Their best-fit 
model requires an H 2 formation rate on dust grains of 3.8 x 
lO^^^cm^ s“' at a temperature of 1000 K, which significantly 
exceeds the standard yalues for H 2 formation on dust grains at 
1000 K. Furthermore, the dust grain FUV extinction cross sec¬ 
tion must be reduced by a factor of ~ 3 while the photoelectric 
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heating rate in the atomic zone requires an increase of a factor of 
~ 3 with respect t o the a priori estimates. 

iHabart et al.l (l2004l) resolved these deviation s from the 
canonical H 2 formation rate coefficient (lJuralll97^ by includ¬ 
ing H 2 formation on very small grains and polycyclic aromatic 
hydrocarbons (PAHs), and adjusting the barriers for H 2 forma¬ 
tion on dust grains. Their model does not fully reproduce the 
observed H 2 formation rates in warm, high-UV environments, 
such as p Oph W and the Orion Bar, where rate coefficients as 
high as 10^'^cm^ s“* have been observed. 

PA Hs make up the low-mass en d of the grain size distri¬ 
bution dWeingartner & Drainell200T^ . The bulk infrared emis- 
sion features of PAHs hav e been detected in interstellar spectra 
dAllamandola et allll989t). but it is still impossible to identify 
individual molecules (IPilleri et alJ l2009h . Although PAHs are 
much smaller than the rest of the dust grains, they account for 
approximatel y half of the total dust grain su rface area available 
for chemistry dWeingartner & Drainel200rB) . In the past decades 
a plethora of studies has been performed to understand how these 
PAHs can catalyze the formation of molecular hydrogen. 

One of the most widely used mechanisms to form H 2 on 
PAHs consists of two different reactions, in which a gas phase 
hydrogen atom forms a covalent bond with a PAH and a second 
step where a different H atom abstra cts a hydrogen atom from 
the PAH to form molecular hydrogen dBauschlichedl998h . Both 
reactions have been studied theoretically as wel l as experimen¬ 
tally w ith sy stems as large as coron ene (C^/iHnl. lMennella et al.l 
d2012h and iThrower et al.l d2Q12h exposed coronene films to 
atomic hydrogen and deuterium to derive cross sections for both 
the addition to and abstraction fror n coronene. Density func- 
tional theory (DFT) calculations by iRauls & Hornekaeif (l2008h 
indicated the presence of barriers when hydro genating a PAH, 
which were later experimentally observed by iBoschman et al.l 
(l2Q12h . These thermal barriers, which have a typical height of 
a few 100 K (a few times 10 meV), significantly restrict the for¬ 
mation of H 2 on PAHs at low temperatures (T < 200K). 

Alternatively, PAH molecules may lose an H 2 molecule af¬ 
ter being excited by a UV photon. This photodesorption mecha- 
nism has been stud i ed both theoretically and in the laborator y 
(lAllain et al.l Il996t iJochims et al.l 1 1994c iLe Page et al.l 120011) . 
These studies find that loss of an H atom is the most dominant 
loss process. However, H 2 loss is only a few times smaller than 
H loss, which could make the former process an important can¬ 
didate for H 2 formation on PAHs, on the condition that H atoms 
also hydrogenate the PAH molecules. It should be noted that the 
relative importance of these processes probably depends on the 
size of the PAH molecule. 

Several theoretical models have been developed to quantify 
the formation of H 2 on PAHs in astrophysical environments. 
These models account for the formati on of H? on PAHs in d iffer- 
ent ways. The Meudon PDR code by iLe Petit et al.l (l2006h con¬ 
siders PAHs as a small-sized extension of the dust grain size 
distribution, where PAH are treated as dust grains with a radius 
between 1 and 50 A. A model bv iLe Page et akl (l2009l) does treat 
the chemical reactions between H atoms and PAHs explicitly. 
However, these reactions are limited to the normal hydrogena¬ 
tion state and a singly dehydrogenated state for both neutral and 
singly ionized molecules. Additionally, it considers a singly su- 
perhydrogenated state for PAH cations. iMontillaud et al.l OQ13h 
presented a model that studies the hydrogenation and charge 
states of PAHs in PDRs using the Meudon PDR code, where 
the possible hydrogenation states range from fully dehydro¬ 
genated up to the singly superhydrogenated state. These hydro¬ 
genation states are available as either a neutral molecule or a 


si ngly ionized cation. Alth oug h the PAH chemistry is different 
in IMontillaud et al.l (1201 3l) and iLe Page et al.l (1200^ . they both 
find that PAHs only exist in either the highest or the lowest pos¬ 
sible hydrogenation state, depending on the column density. In¬ 
termediate PAHs are rarely obseryed. 

In this study we model the chemistry of PAHs in photodis¬ 
sociation regions (PDRs) to inyestigate their ability to form H 2 . 
We combine H 2 formation through LH and ER on dust grains 
and through abstraction and photodesorption on PAHs to estab¬ 
lish total H 2 formation rates at different depths in a PDR. We 
yary the number density of the cloud and the intensity of the 
external radiation field of the PDR. This allows us to study the 
effect of the radiation field on the PAH chemistry. 

We use a time-dependent rate equation model and report the 
results obtained at steady state. Since addition barriers haye only 
been deriyed for coronene, we use this molecule as our generic 
PAH. The transferability of these results is discu ssed in Sect. 14.21 
We als o incorporate the barriers determined in IBoschman et al.l 
(l2012h for the addition and abstraction of H atoms to coronene, 
which prohibits the chemical hydrogenation or abstraction cycle 
for H 2 formation on PAHs at gas temperatures below 200 K. In 
this manner, we find a high-temperature pathway for the forma¬ 
tion of H 2 . 

In Sect. |2] we describe the details of the physics and chem¬ 
istry of our model. The results of the model are shown in Sect. 
|3] The uncertainties of the model are discussed in Sect. 01 where 
we also compare our results and preyious studies. 


2. Model 


We describe a cloud by a one-dimensional PDR model with 
constant density n^. The cloud is irradiated from one side by 
a star with a UV intensity field Go, denoting the strength of 
the radiation field between 6 < h v < 13.6 eV i n Habing units 
(1 Go = 1.6 ■ 10“^ ergcm^^ s~*. lHabinglll968l) . The cloud is 
diyided into equally sized parcels of gas, which are consecu- 
tiyely brought into chemical equilibrium. The photoprocesses 
that are included are photoionization of PAHs and carbon atoms, 
photodissociation of H 2 , and photodesorption of H and H 2 
from PAHs. All photoprocesses are mitigated by dust extinc¬ 
tion at increasing column densities, and for H 2 photodissocia¬ 
tion we also include H 2 self-shielding. Ionization processes are 
coun teracted by electro n recom bination, for which we use rates 
from RuiterkamB_eL^ (l2005l) and the online database KIDA 
(IWakelam et al.l 20 121 http: //kida. obs. u-bordeauxl. £r). 
The formation of H 2 on interstellar dust grains (silicates and 
amorphous carbon) as well as on PAHs is considered in our 
model. 


After a parcel has reached chemical equilibrium, we correct 
the radiation field for dust absorption and H 2 self-shielding for 
the next parcel. The process is then repeated until the entire cloud 
is in chemical equilibrium. Dust grain and coronene abundances 
scale linearly with the number density. The parameters defining 
the cloud are summarized in Table [1] 

The processes inyolying PAHs that we consider are detailed 
in the following subsections. In Sect. l2.1.Tl we describe the pro¬ 
cesses of hydrogenation and abstraction of H atoms on PAHs. 

In Sect. I2.1.2l we discuss the interactions of PAHs with the 
impinging radiation field, where we describe the process of pho¬ 
toionization and electron recombination. We also describe the 
process of photodesorption, where a hydrogen atom or an H 2 
molecule is remoyed from the PAH upon absorption of a photon. 
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Table 1. Model parameters 


Parameter 

Definition 

Value(s) 

Go 

radiation held 

IQiU- 10 ^-^ 

wh 

number density 

10 ^ - lO^cm-3 

^ add 

hydrogenation cross section 

1.10 aVc atom 

^ abs 

abstraction cross section 

0.06 A^ / C atom 

Eadd 

hydrogenation barrier 

see TablelA.il 

Eabs 

abstraction barrier 

lOmeV 

NH,max 

maximum column density 

lO^^cm^^ 

^ph 

photodesorption rate 

see TablelA.2l 


This process is not limited to fully aromatic coronene 
(C 24 H 12 ), but can also occur for coronene in other hydrogena¬ 
tion states, with up to 24 additional hydrogen atoms attached to 
the coronene (C 24 H 36 ). 

A (hydrogenated) coronene molecule can be either neutral or 
positively or negatively charged, which significantly changes the 
chemistry. We do not consider the PAH anions since their contri¬ 
butions are mostly important in strongly shielded environments, 
but concentrate on the barely shielded edges of the PDRs. The 
system of chemical reactions between coronene and hydrogen is 
summarized in Eqs.[T]and|2and is detailed in Table lA.il 


The influence of dust grains on the PDR is described in Sect. 
12.21 and we treat the photodissociation of H 2 in Sect. 12.31 In Sect. 
I2.4l we discuss the model calculations themselves. 


2.1. PAHs 

The abundance of PAHs in the interstellar mediurn (ISM) i s very 
uncertain, and we refer to IWeingartner & Drainel (l2001ah for a 
comparison between different values in several environments. 
We assume a total carbon abundance of 2.5 ■ 10 of which a 
third i s loc ked in PAHs. This is similar to studies bv iJoblin et al.l 
(Il992h and lLi & Drainel (1200 ih . who find the fraction of carbon 
atoms locked in PAHs relative to hydrogen to be 6 x 10“^. 

We define PAHs as being smaller than 100 A, and we 
use coronene (C 24 H 12 ) as our prototypical PAH of the ISM. 
This is motivated by the fact that coronene is a well-studied 
molecule, both theoretically and experimentally. The barriers for 
single H addition and abstraction, which are instrumental for 
this study, have been determined only recently. In this sense, 
we use coronene as a molecule representative of all the PAHs 
in the PDRs studied in this work. It should be noted that this 
is a rough assumption, since it is widely believed that small 
PAHs (Nc < 40) cannot survive in strong radiation fields, 
which would exclu de the existence of the coronene molecule in 
PDR environments (lAllain et alJll996t) . an d that compact PAHs 
are more stable t han extended molecules (iJochims et al.l 1 19^ 
lEkern et alJl997h . However, other studies show that small PAHs 
can s urvive strong radiation f ields, but in a very dehydrogenated 
state (iMontillaud et al.ll2013h . Therefore a more realistic PDR 
model should consider a distribution of PAHs with different 
sizes, where larger PAHs with fewer H-bearing outer edge car¬ 
bon atoms would be present. This implies that these PAHs would 
be slightly less efficient in forming H 2 , with an efficiency de¬ 
pending on the number of inner carbon atoms compared to outer 
edge carbon atoms. Although this would not affect the results 
drawn in this work, in future work a size distribution of PAHs 
would further increase the reliability of our model. In Sect. 14.21 
we discuss the expected results when coronene is replaced with 
such a size distribution. 


2.1.1. Hydrogenation and abstraction 

Hydrogen addition and subsequent abstraction constitute a two- 
step catalytic cycle for the formation of H 2 . Eirst, a coronene 
molecule reacts with a hydrogen atom to form a singly super- 
hydrogenated coronene molecule (C 24 H 13 ). In the second reac¬ 
tion, another hydrogen atom reacts with this singly superhydro- 
genated coronene mol ecule, yielding mole cular hydrogen and 
the coronene molecule (lBauschlicheiill99^ . 


c24h°/+ + h ^ c24h;;^; (1) 

c24h°^; + h ^ c24h;;/+ + h2. (2) 

These reactions take place at a rate r, which is calculated as 
r(n) = n (H) n (C 24 H“^+) cry* exp j ■ (3) 

Here, n (H) and n (C 24 Hn^^) represent the number densities 
of atomic hydrogen and coronene molecules or cations contain¬ 
ing n H atoms, respectively, cr represents the cross section for 
the specific reaction and y* is the average thermal velocity of 
atomic hydrogen. We use a cross section of 1.1 A^ per reactive 
carbon atom for hyd rogenation and 0.06 A ^ per reactive carbon 
atom for abstraction (iMennella et al.ll201^ . 

Since coronene is much more massive than hydrogen, its ve¬ 
locity is neglegible compared to that of the gas phase hydrogen 
atoms, y* is therefore the thermal velocity of hydrogen atoms, 
which is calculated as 



(4) 


where ke is the Boltzmann constant, T is the gas temperature 
and oth is the mass of the hydrogen atom. 

The hydrogenation and abstraction of coronene is subject to 
barriers (E^), depending on hydrogenation and charge state. The 
presence of these b a rriers was predicted for neutral coronene by 
iRauls & Hornekae^ (!2008h and first experimen tally determined 
for coronene cations by lBoschman et akl (l2012l) . The latter study 
shows an alternating behayiour for the existence of barriers, 
meaning that the first hydrogenation of the coronene cation 
is barrierless. The second hydrogenation then has a barrier, 
while the third does not, etc. More r ecent experiments in con- 
juntion with theoretical calculations (iBoschman et al.llin DreD.h 
show that this alternating behayiour continues until the coronene 
cation is fully hydrogenated. The hydrogenation barriers are 
summarized in Table I A. fl 

Abstraction barriers are only present for abstraction from 
those hydrogenation states for which adding a hydrogen atom 
would also haye a barrier. In this manner, a distinction is made 
between chemically stable and chemically unstable hydrogena¬ 
tion states, where the unstable ones haye a radical atom. The 
barrier again st abstraction is set at 10 meV, equal to the yalue 
calculated by lRauls & Hornekaed (l2008l) for abstracting an outer 
edge H atom from the doubly hydrogenated coronene molecule 
(C 24 H 14 ). The influence of the uncertainties of these barrier 
heights on the obtained results is discussed in Sect. |4] 
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With the mechanism for H 2 formation catalyzed by PAHs as 
described in Eqs. ([TJ and (| 2 ]l, the total H 2 formation rate due to 
abstraction from PAHs is calculated as 

36 

/?/(PAHabs) = ^r(n), (5) 

n=0 

where r(n) is the H 2 formation rate as shown in Eq. ([3]). 


with T the temperature, s (e) the electro n sticking coefficient an d 
We the number density of free e lectrons (iRuiterkamn et al.l2005h . 
Erom iRuiterkamp et al.l (|2005|) we adopt an electron sticking co¬ 
efficient of s (e) = 10 “^ for coronene. 

We only use carbon atoms as a source of free electrons, be¬ 
cause other metals with an ionization energy comparable to car¬ 
bon have abundances that are at least one order of magnitude 
lower. Since the ionization balance is driven by radiation, other 
metals are not considered. 


2.1.2. Photodesorption 

In a PDR, PAHs experience a strong interaction with the pho¬ 
tons from the impinging radiation field, which can either ionize 
the PAH or excite i t into a higher state. We use Eq. dSll from 
iLe Page et al.l (l2001l) to calculate the total photoionization rate: 

^13.6 eV 

^ion= YUE)cTvw{E)N{E)dE, (6) 

Jip 


with IP the ionization pot ential of the molecu le and Fion the 
ionization yield taken from iLe P age et a ll (12001 ). cruv(E) is the 
photoabsorption cross section JMalloci etak 2007h and N{E) is 
the intensity o f the Habing interstellar radiation field in photons 
cm^^s^'eV ' (lDrainelll97^ . 

If the molecule is not ionized but excited, it can dissipate the 
excitation energy by ei ther IR emission, loss of an H atom, or 
loss of an H 2 molecule (lAllain et al.ll996ll . To calculate the rates 
for these processes, we use the RRKM method as described in 
iLe Page et al.l (l2Q0ll) . which gives the rate of bond dissociation 
as a function of excitation energy, bond energy and vibrational 
degrees of freedom in the molecule. The competition between 
processes is accounted for by dividing the rate of one process by 
the sum of the rates of all possible processes. In this way, the 
branching ratio for each dissociation pathway upon excitation 
is obtained. Multiplying this dissociation probab ility with the 
corone ne photoabsorption cross section cruv(E) of iMalloci et alJ 
(l2007h and integrating over the Habing interstellar radiation field 
N(E) below the Lyman limit gives the total photodissociation 
rate per Gq, as shown in Eqs. O and ([ 8 ]): 


/' 


^13.6 eV 

/:ph(H)= I (l-Fion)- 


/:diss(H) 


^diss(H) + + /^diss(H2) 

X o-uv(E)N(E)dE\ (7) 

Jo "" koissCW + k'^^ + koiM) 

X o-uv(E)N(E)dE, (8) 


where is the rate of IR stabilizatiom _ 

We us e bindi ng energies from iLe Page et al.l (1200 Ih and 
iParis et aTI (l2014l) to calculate the photodissociation rates for 
both the neutral coronene and the coronene cation in their de¬ 
hydrogenated, superhydrogenated, and normally hydrogenated 
states. These photodissociation rates cover both loss of a single 
H atom and the loss of H 2 and are summarized in Table lA.2l 
In addition to PAH ionization, we also include the photoion¬ 
ization of unbound carbon atoms as a source of free electrons. 
Recombination reactions are included in our chemical network, 
which ensures an ionization balance. The rates for recombina¬ 
tion reactions on PAH cations are given by 


k 


rec 


= 3 X 10 ^^ 



s(e) Me, 


(9) 


2.2. H 2 formation on dust grains 


Dust grains have a significant influence on the chemistry taking 
place in the PDR because of their catalyzing properties, which 
strongly depend on the chemical composition of the dust grain, 
for example, silicates or amorphous carbon. 

In our model, we limit the dust grain chemistry to the catal¬ 
ysis of H 2 formation, and for this we assume the first H atom to 
be either physisorbed or chemisorbed. There are then two mech¬ 
anisms for this H atom to react with another H atom and to form 
H 2 . In the first process, an H atom from the gas phase interacts 
with an adsorbed H atom (Eley-Rideal). In the second process, a 
different physisorbed H atom travelling around the surface inter¬ 
acts with the adsorbed H atom (Langmuir-Hinshelwood). 

Eor both mechanisms to occur, a hydrogen atom must first 
adsorb onto the surface of the dust grain and become either 
chemisorbed or physisorbed. This adsorption process is incor¬ 
porated into rat e equations with a tempera ture-dependent stick¬ 
ing coefficient (iBurke & Hollenbachlll983ll . Once an H atom is 
adsorbed onto the surface, it has a certain probability to form 
H 2 . This probability e strongly depends on the binding energies 
of the adsorbed atoms and the gas and dust temperatures. The 
resulting H 2 formation rate equation is given by 


^•giam = n (grain) n (H) Uthcre^. 


( 10 ) 


Here n (grain) is the number density of dust grains, n (H) is 
the number density of atomic hydrogen, and t^th is the aver¬ 
age thermal velocity in the gas phase, cr is the average cross 
secti on of a dust grain, and values for n (grain) cr are taken 
from ICazaux & Snaansl (l2009h . Eor the values of formation ef- 
ficiency e and sticking c oefficient S , we use the expressions of 
ICazaux & Sp aans ( 20091) and the dust g rain size distribution de- 
riyed by lWeingartner & Drainel(l2001a ). The sticking coefficient 
is calculated as 


S (Tgas, Tdust) — 


1 -H 0.4 ■ 


- gas 


Tdu 


100 


0.5 


-I- 0.2' 


100 


-H 0.08 ■ 


T \ 

^ gas I 

100 / 


23 


( 11 ) 


which giyes a numerical yalue for S typically between 0.4 
and 0.9. The formation efficiency strongly depends on the sur¬ 
face composition and has a typical yalue between 0.4 and 1 
for amorphous carbon; i t lies between 0.01 and 1 for silicates 
(ICazaux & SDaansll 2009 ). _ 

ICazaux & Snaansl (l2009h considered H 2 formation on the 
surfaces of silicates, amorphous carbon and PAHs. We ayoid a 
redundancy by only taking into account H 2 formation on sili¬ 
cates and amorphous carbon and defining this as H 2 formation 
on dust. The H 2 formation on dust is calculated by using the sur¬ 
face area of silicates ( 1 ^ = 10 ^^'cm^^) and amorphous carbon 
/ Vf; — YJ X 10 “^^cm“^), so as to not count the contribution of 
PAHs twice. 
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2.3. H 2 photodissociation 

The destruction mechanism for molecular hydrogen in the pres¬ 
ence of UV photons is dissociation subsequent to the photo¬ 
excitation into a higher electronic state. Approximately 10% of 
these excited H 2 molecules undergo fluorescence into the vibra¬ 
tional continuum of the electronic ground state, giving rise to 
dissociati on of the molecule. Th e H? p hoto dissociation rate is 
given by iBlack & van Dishoecld (119871) and iDraine & Bertoldil 
(119961) : 


kdiss = 3.4- 10-'7(A(H2))Goexp[-2.6Av]s-'. (12) 


In this equation the exponential factor expresses the dust ex¬ 
tinction. The factor /(A^(H 2 )) accounts for H 2 self-shielding 
against photodissociation. This self-shielding factor /(A(H 2 )) 
has been computed as a function of H 2 column density by 
iDraine & Bertoldil d 19961) . 


/(A^(H2)) = 


0.965 
(1 + x/bsf 


0.035 

Vl + y 

X exp [-8.5 ■ 10“"^ Vl + , 


(13) 


where x = A(H 2)/5 ■ 10 '^ cm“^ is the H 2 column density and bs 
is the line width of H 2 absorption lines expressed in 10 ^cm s ', 
which was set to unity for our calculations. 


2.4. Computations 


In our calculations, we use a one-dimensional plane-parallel slab 
of gas that is irradiated from one side. This slab is divided into 
10 ^ equally sized cells of gas, which all have the same num¬ 
ber density and abundances. Initially, there is no molecular hy¬ 
drogen present, and the only PAH is the normally hydrogenated 
coronene molecule (C 24 H 12 ). 

Every single reaction in the list of reactions gives rise to a 
rate equation, and from this the rate for every reaction is calcu¬ 
lated at f = 0. We perform a numerical Euler integration until 
the system reaches convergence, which is defined as the moment 
when the fractional difference between two time steps is smaller 
than 10“^ for every chemical species. Eor the species with abun¬ 
dances higher than 10 this condition is met in less than 10 ^ 
years. Eor species with lower abundances it is increasingly diffi¬ 
cult to converge, and timescales ranging between 10 ^ and 10 * 
years are required to reach equilibrium for the entire system. 
Time-dependent environmental effects may play a role, but this 
will be a topic of future study. After convergence is reached in 
a parcel of gas, the calculation starts for the next parcel of gas, 
including a modification of the radiation field to account for dust 
extinction and H 2 self-shielding. 

Temperatures in a PDR are strongly dependent on the 
external radiation field through various mechanisms. In re¬ 
gions where dust grains are exposed to UV radiation, the 
gas is indirectly heated via photoionization of dust grains 
and PAHs. At low densities (nu ^ 3 ■ let'em *), the ther¬ 
mal coupling between dust grains and the gas phase does 
not occur, which leads to gas temperatures significantly 
different from the dust temperatures. At higher densities 
(hh ^ 3 ■ lO^cm^*) thermal coupling between gas and dust 
grains has to be taken into account. We extract gas and dust 
pain t emperatures from models of PDRs (iMeiierink & Snaani 
20051) and refined by iHocuk & Snaansl ( 2011 ). available at 


http://www.mpe.mpg.de/~seyit/Tgas/). These models 
derive gas and dust grain temperatures from the different heating 
and cooling processes, leading to temperatures for a given 
combination of radiation field, number density, and column 
density. We use the outcomes of these models in tabular form, 
from which we extract the temperatures for our simulations. 
Temperature fluctua tions of du s t grai ns are not considered in 
this model, because iBron et aTI (12014h showed that their effect 
on the H 2 formation rate is limited. 

The external radiation field has a significant influence on the 
structure of the cloud through photoreactions and heating of the 
gas. At high total column densities (Nh ~ 10^' cm“^) and at solar 
metallicity, the intensity of the radiation field diminishes due to 
dust absorption and scattering. The radiation field is therefore a 
function of column density, as shown in Eq. [14] 

Go(Av) = Go(0)exp[-/? Ay]. (14) 

In this expression, fi is the process-dependent extinction fac¬ 
tor, which has a typical value of 2.5, and Ay is the visual ex- 
tinction calculated b y Ay = (V(H) -H 2 V(H 2 ))/ 2.2 ■ lO^'cm-^ 
(iGuver & Ozelll2()()^ . 

3. Results 

Erom our model we obtain number densities for all coronene 
charge and hydrogenation states, as well as the H 2 formation 
rate associated with each of these species. We also obtain the 
H 2 formation rates on dust grains and the number densities of 
atomic and molecular hydrogen. Together with the gas and dust 
grain temperatures, this gives us a clear view of the chemistry 
and physics inside the PDR. We present the spatial distribution 
of coronene hydrogenation states in Sect. 13.11 In Sect. 13.21 we 
show the H 2 formation rates throughout the PDR, and in Sect. 
I3.3l we show how the formation of H 2 on coronene impacts the 
PDR structure. 

To study the hydrogen-coronene chemistry, the spatial dis¬ 
tribution of coronene hydrogenation states is tracked. A region 
with low hydrogenation states (C^ and C 24 H''’^ dominate) in¬ 
dicates that photodesorption occurs at a rate significantly higher 
than that of hydrogen addition. Conversely, in a region with nor¬ 
mal and high hydrogenation states (C 24 Hj 2 ''' and higher domi¬ 
nate), hydrogenation dominates H-loss processes. Since H 2 for¬ 
mation is possible through both photodissociation of coronene 
and chemical abstraction on coronene, the spatial distribution of 
the hydrogenation states of coronene can tell us which of these 
processes is more important for H 2 formation. Furthermore, in 
Sect. l3.2l we compare the formation of H 2 on coronene with the 
formation on dust grains, and report our results. 

3.1. Spatiai distribution of coronene 

The spatial distribution of coronene states is due to the tempera¬ 
ture and radiation profiles of the cloud. At the boundaries of the 
cloud, the radiation field is very strong and the gas temperature is 
high (Tgas > 2QQK), as is shown in Fig.[T] At these high gas tem¬ 
peratures coronene hydrogenation is very efficient because the H 
atoms can then easily overcome the barriers against hydrogena¬ 
tion. However, this is easily offset by the intense radiation field 
that completely dehydrogenates coronene down to the carbon 
skeleton. Deeper inside the cloud the radiation field is decreased, 
causing a decrease in the photoionization rate and a subsequent 
gradual transition from fully dehydrogenated coronene cations 
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Fig. 1. Distribution of dust and gas temperatures throughout a gas cloud 
with Go = 10^'^ and uh = 10^'^cm^^. The dotted lines show the distri¬ 
bution of H and H 2 . 


to their neutral counterparts is expected. At even higher column 
densities, the radiation field is too weak for photodissociation to 
dominate. As a result, a transition to the neutral, normally hy¬ 
drogenated coronene molecule should occur. 

A typical spatial distribution of hydrogenated coronene is 
presented in Fig.|2 In the outskirts of the cloud coronene is fully 
dehydrogenated and exists as cations. Deeper inside the cloud, 
the dominant form of coronene is the fully dehydrogenated neu¬ 
tral molecule, which is subsequently overturned by the normally 
hydrogenated neutral coronene molecule. 

Clearly absent in our results are the highly superhydro- 
genated species, although coronene with a few extra hydrogen 
atoms can be abundant at high extinctions. This can be explained 
by the cori'elation between the UV intensity and the gas tem¬ 
perature. A high gas temperature is necessary to overcome the 
energy barriers in the hydrogenation process. These gas tem¬ 
peratures are only reached in regions where the impinging ra¬ 
diation held is strong enough to heat the gas to such tempera¬ 
tures. At the same time, this radiation held is strong enough to 
photodissociate all the coronene molecules. This is consistent 
with p revious PDR models, such as those of iMontillaud et akl 
( I2OI3I) . who modelled the north-west PDR of NGC 7023. They 
reported that at every depth there is one dominant hydrogenation 
state, which is either the one with the most H atoms removed 
or the highest possible hydrogenation state. Their hydrogenation 
state s range from full dehydr ogenation to one superhydrogena¬ 
tion. IMontillaud et al.l (1201 3li concluded that coronene is fully 
dehydrogenated throughout the entire cloud (Ay < 3 mag). This 
is consistent with our findings, where we report that coronene 
is dehydrogenated at the edge of the cloud, but normally hydro¬ 
genated at extinctions of more than a few mag. 

The results shown here are also consistent with a recent study 
that interpreted Sp itzer and ISOCAM observations of NGC 7023 
( Joblin et al.ll2010l) . Using the method developed by Pilleri_et^ 

( 20121) to identify PAH cations and neutral PAHs, iJoblin et^ 
(I 2 OIOI) showed that there is a sharp transition between PAH 
cations at the edge of the PDR and neutral PAH molecules deeper 
into the cloud. 

The radiation field has a very clear impact on the coronene 
distribution through two different mechanisms, as shown in Fig. 
|2] First, an increased radiation field leads to a higher H 2 pho¬ 


Fig. 2. Distribution of coronene charge and hydrogenation state for three 
different values of Go and rtu- Coronene cations are plotted with dashed 
lines, and neutral coronene molecules have solid lines. Only the most 
abundant hydrogenation states are shown. The top graph has Go = 10“'^ 
and wh = 10^'^cm^^, the middle graph Go = 10^° and uh = 10^'^cm^^ 
, and the bottom graph corresponds to Go = 10^ ^ and nu = lO'^'^cm^^. 
The vertical line marks the location of the H/H 2 transition, where «(H) = 
2 • n(H2). 

todissociation rate, pushing the H/H 2 front to higher column den¬ 
sities. Second, a higher Go leads to a higher rate of photoioniza¬ 
tion of PAHs, changing the transition from cations to neutrals to 
higher column densities. Furthermore, a stronger radiation field 
drives the transition from dehydrogenated to normally hydro¬ 
genated PAHs to higher column densities. 

An increase in the number density of the cloud has the op¬ 
posite effect of increasing the radiation field. Hence, the hydro¬ 
genated/dehydrogenated front is shifted towards lower column 
densities, as shown in Fig.|2l 

3.2. H 2 formation rates 

The H 2 formation rates on coronene and dust grains are strongly 
dependent on the radiation field and the dust temperature, re¬ 
spectively. For H 2 formation on dust, the dust must be cold 
enough to prevent physisorbed H atoms from desorbing into the 
gas phase. When Tjust rises above 15-20 K, H 2 formation is only 
possible with chemisorbed H atoms and is less efficient. 

For coronene, however, there are two possible pathways for 
H 2 formation, as described in Sects. IZLTI and 17. 1.21 The dom¬ 
inant H 2 formation process is the photodesorption of H 2 from 
coronene with at least two hydrogen atoms. Although the oc¬ 
curence of these molecules is low in PDRs, this process occurs 
very fast because of the intense radiation. This results in a fast 
catalytic cycle where the bare carbon skeleton of the coronene 
molecule consecutively adsorbs two H atoms, which are sub¬ 
sequently photodesorbed as an H 2 molecule. The H 2 forma¬ 
tion rates for this mechanism can be as high as a few times 
10~*^cm ^s~*, which i s simil ar to H 2 formation rates in PDRs 
found bv iHabart et al.l (l2004l) . 

The formation of H 2 through chemical abstraction on 
coronene is not a dominant process in PDR environments, al¬ 
though the temperatures are high enough to easily cross the ther¬ 
mal batTiers associated with this process. For the latter process 
hydrogen atoms must be adsorbed to the coronene molecule, but 
the high-radiation field photodesorbs almost all H atoms from 


Article number, page 6 ofM 












































Boschman et al.: Hydrogenated PAHs in PDRs 


I — Cations - ■ Neutrals - - Dust| 



0.0 0.2 0.4 0.6 0.8 l.i° 

JV„ (10^^ cm"^) 


Fig. 3. H 2 formation rate coefficients on cationic and neutral coronene 
and dust (solid, dashed-dotted, and dashed lines, left y-axis) and dis- 
trihutions of H and H 2 (dotted lines, right y-axis). The top graph has 
Go = 10^'^ and nn = lO^^cm^^, the middle graph Go = 10'^® and 
«H = 10^ ^cm“^ , and the bottom graph corresponds to Go = 10^'^ and 
«H = lO'^'^cm^^. The vertical line marks the location of the H/H 2 transi¬ 
tion, where )i(H) = 2 • n(H 2 ). 


the coronene molecule. This leaves most coronene molecules 
without H atoms, and a small fraction has only one or two hy¬ 
drogen atoms, rendering the chemical abstraction process negli¬ 
gible. 

To compare H 2 formation on coronene with H 2 formation 
on dust grains, we present the H 2 formation rate coefficients as 
a function of column density in Figs. [3] and |4] In the same hg- 
ures the number densities of atomic and molecular hydrogen are 
reported. To study the effect of the radiation held and the total 
hydrogen number density, the simulations presented in the three 
panels of Fig. [^differ in radiation held or number density hh- 

At low column densities, the gas and dust are heated by the 
impinging radiation held, making coronene the dominant source 
of H 2 formation. At increasing column densities, the tempera¬ 
tures start to decline, and dust grains take over as the dominant 
H 2 formation agent. 

Increasing the radiation held heats the gas and the dust 
grains, making the transition from PAHs to dust grains occur at 
higher column densities, which is shown in the middle panel in 
FigE] Increasing the density has the opposite effect of increasing 
the radiation held, causing a decrease of the column density for 
the PAH/dust grain transition, as illustrated in the bottom panel 
of Fig. a 

3.3. Impact of H 2 formation on coronene on the PDR 
structure 

H 2 formation on coronene can signihcantly change the structure 
of a cloud. This inhuence manifests itself in the column density 
of the H/H 2 transition in the cloud. We consider the H/H 2 tran¬ 
sition to occur at a column density where 2n(H2) = n(H). This 
column density is lowered by the presence of PAHs in the PDR. 

An example is presented in Fig. 01 where we show the num¬ 
ber density prohles of atomic and molecular hydrogen. In the 
top panel of Fig. 01 H 2 can be formed by reactions on coronene 
and dust grains, and in the bottom panel only H 2 formation on 
dust grains is considered. Comparing these hgures, we see that 
the column density of the H/H 2 transition is lowered by a factor 


of 4 due to the formation of H 2 on coronene. The influence of 
coronene can be quantihed by the ratio of the column densities 
of the transition with and without considering the formation of 
H 2 on coronene. We dehne the influence parameter /pah as 


_ iv (//);,,„, 

^PAH - i ” 


dust + PAHs 


dust 




(15) 


which ranges between 0 and 1. In the case of a value of 0, 
coronene has no influence on the column density of the H/H 2 
transition. When /pah = 1, coronene shifts the H/H 2 transition 
completely to the edge of the cloud. The highest value of /pah 
present in our results is approximately 0.8. At this value, the rate 
of H 2 formation on coronene is 200 times higher than H 2 forma¬ 
tion on dust grains at the column density of the H/H 2 transition, 
while these rates are approximately equal for the lowest values 
of the influence parameter. This influence parameter is shown as 
a function of cloud number density and external radiation held 
in Fig.|5] The inhuence parameter increases with both increasing 
number density and external radiation held. 

An increased number density increases the rate of reactions 
between coronene and H atoms, resulting in a higher abundance 
of PAHs with one or more hydrogen atoms. At a similar radiation 
held, this yields a higher H 2 formation rate, pushing the H/H 2 
transition to lower column densities. A similar argument applies 
to the dependence on the radiation held. With an increasing radi¬ 
ation held, all photodissociation rates increase, including those 
of the H 2 forming processes. However, the effect is more subtle 
because a stronger radiation held changes the population dis¬ 
tribution of the hydrogenation states towards a lower fraction of 
coronene with two or more H atoms. This makes fewer coronene 
molecules eligible for the release of H 2 upon photoabsorption. 
This effect leads to an optimal radiation intensity for every num¬ 
ber density for maximal coronene inhuence. If the radiation held 
has a higher intensity, the inhuence of coronene on the position 
of the H/H 2 front decreases again, as is shown in Fig. |5] 

Strong radiation helds also correlate with high dust grain 
temperatures, which inhibit H 2 formation. Figure|6]shows the H 2 
formation rate of dust grains (solid lines) and coronene (dashed 
lines) as a function of gas temperature for different values of Go 
and hh- For each datapoint, the color of the marker indicates the 
dust tempe rature. Additionally, observational H 2 formation rates 
derived by [Habart et al.l (l2004l) for different PDRs are reported 
in this hgure. It is clear that H 2 formation on PAHs is necessary 
to reproduce the observed H 2 formation rates. 

From this plot it is also clear that dust grains and PAHs con¬ 
tribute to H 2 formation in different temperature regimes. Since 
the temperature decreases with increasing column density, PAHs 
have their strongest H 2 formation activity at the edges of PDRs, 
while dust grains dominate this process deeper inside the cloud. 
PAHs change the cloud structure the most when the conversion 
from H to H 2 occurs at lower extinction than the one at which 
dust grains start to dominate H 2 formation. We can see from Fig. 
|6]that this is the case when the transition occurs in a region where 
Tgas > 200 K. 

Combining high number densities and radiation helds, H 2 
formation on PAHs builds up a sufficiently high H 2 column den¬ 
sity at a relatively low total column density, such that H 2 is al¬ 
ready self-shielded against photodissociation. At these low total 
column densities, dust temperatures are still too high for the for¬ 
mation of H 2 on dust grains to have any effect (Tdust > 40K). 
The column density of the H/H 2 transition is in that case entirely 
due to H 2 formation on PAHs. 
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Fig. 4. Comparison between two identical clouds, with (top panel) and 
without (bottom panel) H 2 formation on coronene. H 2 formation rate 
coefficients on ionized (solid) and neutral (dashed-dotted) coronene and 
dust grains (dashed lines) are represented on the left y-axis. Number 
densities of H and H 2 (dotted lines) are represented on the right y-axis. 
The thick vertical green line indicates the location where n(H) = 2 • 
n(H 2 ). Go = lO'^'^ and nn = lO^cm^^. 
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Fig. 5. PAH influence parameter plotted as a function of number density 
and external radiation field. The influence of PAHs on the cloud struc¬ 
ture increases with increasing number density and radiation field. White 
regions indicate that no data are available. 

4. Discussion 

4. 1. Impact of uncertainties 

We find that PAHs, with coronene as a prototypical PAH, affect 
the structure of a PDR. The formation of H 2 on coronene is pri¬ 
marily important on the outskirts of the PDR. In these regions 
the temperature is sufficiently high that the exact height of the 
barriers for addition and abstraction does not matter. Deeper in¬ 
side the cloud the barrier height does matter, but the transition to 
molecular hydrogen is then already complete, so H 2 formation 
on PAHs is not significant in these regions. 

Furthermore, most of the H 2 formation on PAHs takes place 
via photodesorption. These processes occur at low column den¬ 
sities, where the radiation field suffers almost no extinction. At 
these column densities, photodesorption processes dominate and 


Fig. 6. Hj formation rates on for different values of Go and «h as a 
function of Tgas- Solid lines indicate H 2 formation by PAHs and dashed 
lines represent dust grain H 2 formation. For each datapoint, the color of 
the marker corresponds to the dust temperature of that datapoint. Both 
the gas and dust tempe ratures have heen derived from the models in 
iHocuk & Snaansl j201l|). The hlack markers indicate the H 2 formation 
rates reported hv iHabart et all | |2004|) for different PDRs. These forma¬ 
tion rates have been estimated using H 2 line intensity ratios from the H 2 
fluorescent lines. 


have a strong influence on the depth of the H/H 2 transition. To 
understand the dependence of the H/H 2 transition on these pho¬ 
torates, we vary the binding energies of either H atoms or H 2 
molecules with 0.2 eV ceteris paribus and recalculate the pho¬ 
torates. With these photorates the simulations are re-evaluated 
and show only a weak influence on the position of the H/H 2 
front. 

Increasing the H 2 or decreasing the H photodissociation en¬ 
ergy (Eq in Table lA.2l i with 0.2 eV results in H 2 formation rates 
reduced by a factor 3. Conversely, the H 2 formation rates locally 
increase by a factor 3 with a similar decrease in the photodissoci¬ 
ation energy for H 2 loss, or an increase in the photodissociation 
energy for H loss. However, this does not significantly change 
the position of the H/H 2 front. Furthermore, we see in our simu¬ 
lations that in the clouds where PAHs have a strong influence on 
the cloud structure, their H 2 formation rate is typically one order 
of magnitude higher than the dust grain H 2 formation rate. 

Another parameter that is subject of discussion is the cross 
sections used for the hydrogenation and abstraction reactions. 
The reaction rate scales linearly with the cross section, and the 
values used in this paper are an average value derived from ir ¬ 
radiating coronene fil ms with H atoms (iMennella et al.llT012l) . 
iMennella et al.l (l2012h argued that because coronene molecules 
only weakly interact, and since the H atoms interact directly with 
the coronene molecules, without diffusion or absorption, their 
values are representative for the gas phase molecules as well. In 
our model abstraction plays a minor role in H 2 formation, but this 
changes for larger PAHs (Sect. l4~2]) . In the cases where abstrac¬ 
tion is the dominant mechanism for H 2 formation, the formation 
rate linearly depends on both the abstraction cross section and 
the number of H atoms adsorbed onto PAHs. The latter in turn 
depends on the cross section for the addition of an H atom. We 
therefore expect that the rate of H 2 formation through abstraction 
linearly depends on both the cross section for addition and ab¬ 
straction. To make more reliable H 2 formation models for larger 
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PAHs than this proof-of-concept approach, both these cross sec¬ 
tions have to be established for larger PAHs. 

4.2. Expectations for a population of large PAHs 

We have used coronene as the only active PAH, where in reality 
there are many more PAHs of different sizes present in the ISM. 
For PDR environments, PAHs with more than 40 carbon atoms 
and a compact structure a re favoured because of their photosta¬ 
bility. lAllain et’aP (1 19961) showed that both H loss and H 2 loss 
significantly decrease with increasing number of carbon atoms. 
This makes H 2 photodesorption a less likely process for H 2 for¬ 
mation, but the decreasing rate of H loss will leave more H atoms 
on the PAH molecule, making H 2 formation through abstrac¬ 
tion a more dominant process. This also changes the distribution 
of hydrogenation states throughout the PDR and increases the 
abundance of higher hydrogenation states. These changes in both 
the dominant H 2 formation mechanism on PAHs and the hydro¬ 
genation state distribution are observed in a test run of our model 
where coronene is substituted by a generic 50C PAH. Since the 
chemistry between PAHs and hydrogen atoms is less well under¬ 
stood for larger PAHs, we did not include these molecules in our 
simulations. 

Larger molecules like circumcircumcoronene (C 96 H 24 ) are 
found to b e in the highest hydroge nation state throughout the en¬ 
tire cloud dMontillaud et al.ir2013h . However, the photo-stability 
of hydrogenated PAHs is quite unknown, and preliminary exper¬ 
imental results suggest that the addition of hydro gen on PAHs 
changes their photostability (IReitsma et al.ll20f^ . This implies 
that the dissociation rates of hydrogenated coronene molecules 
could be lower than previously assumed. As a result, the forma¬ 
tion of H 2 through photodesorption will decrease as well. The 
photodesorption H 2 formation rates presented in this paper are 
therefore an upper limit to possible H 2 formation rates through 
this mechanism. Simultaneously, H 2 formation through abstrac¬ 
tion will become more dominant on larger PAHs, since the de¬ 
creased photodesorption rate leaves more H atoms on the PAH 
molecule. However, to understand to which order abstraction can 
replace photodesorption as an H 2 forming mechanism, a more 
comprehensive understanding of PAH (photo)chemistry is nec¬ 
essary. 

Furthermore, there are other processes for H 2 formation on 
PAHs, such as dissociative recombina tion between a superh ydro- 
genated PAH cation and an electron (iLe Page et al.ll20()^ . This 
process is found to be effective for superhydrogenated PAHs 
with up to 200 carbon atoms. In the environments we study the 
rate for this process is at least an order of magnitude lower than 
the dominant H 2 formation processes. To conclude about H 2 for¬ 
mation on PAHs, it is necessary to fully understand every con¬ 
tributing process, since different processes appear to be the dom¬ 
inating mechanism for different PAH sizes and hydrogenation 
states. 


5. Conclusions 

The model we presented clearly shows that PAHs play a vi¬ 
tal part in the location of the H/H 2 transition in PDRs, espe¬ 
cially at high gas temperatures (T > 200 K). For high UV fields 
(Go > 10^), PAH (photo)chemistry should therefore be taken 
into account when modelling the structure of a PDR. 

We presented the first chemical model that accounts for more 
than one superhydrogenation of a PAH molecule and that incor¬ 
porates thermal barriers into the PAH chemistry taking place. 


This model also considers photodesorption as a mechanism for 
H 2 formation on PAHs. We showed that the barriers introduce 
a strong temperature dependence for the H 2 formation route on 
PAHs, and that this formation pathway is the high-temperature 
complement to the classical H 2 formation mechanism on dust 
grains. Further research into PAHs should be directed at a bet¬ 
ter understanding of dehydrogenated PAHs and the chemical and 
photoprocesses that are associated with them. Ideally, this is also 
directed towards larger PAHs, which are more representative of 
the ISM. 

We conclude that photodesorption of H 2 from PAHs can be 
an important pathway for the formation of H 2 in PDR environ¬ 
ments. The barriers present in PAH hydrogenation limit the H 2 
formation process on PAHs to high gas temperatures. These bar¬ 
riers therefore ensure that H 2 formation on PAHs is the high- 
temperature complement to dust grain H 2 formation. 
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Table A.l. Chemical reactions with different hydrogenation and charge states of coronene, as shown in the first column. The second and third 
column show the barrier and number of reactive carbon atoms for addition of a hydrogen atom. The fourth and fifth column show these properties 
for the process of abstraction, which leads to H 2 formation. 



Hydrogenation 


Abstraction 


Reactant 

Barrier / meV (K) Reactive sites 

Barrier / meV (K) Reactive sites 

C24H0 

- 

12 

- 

- 

C24H1 

- 

11 

10(116) 

1 

C24H2 

- 

10 

10(116) 

2 

C24H3 

- 

9 

10(116) 

3 

C24H4 

- 

8 

10(116) 

4 

C24H5 

- 

7 

10(116) 

5 

C24H6 

- 

6 

10(116) 

6 

C24H7 

- 

5 

10(116) 

7 

C24H8 

- 

4 

10(116) 

8 

C24H9 

- 

3 

10(116) 

9 

C24H10 

- 

2 

10(116) 

10 

C24H11 

- 

1 

10(116) 

11 

C24H12 

28 (324) 

12 

10(116) 

12 

C24H13 

- 

1 

- 

1 

C24H14 

46 (533) 

2 

10(116) 

2 

C24H15 

- 

1 

- 

1 

C24H16 

64 (742) 

1 

10(116) 

4 

C24H17 

- 

2 

- 

1 

C24H18 

35 (406) 

1 

10(116) 

6 

C24H19 

- 

1 

- 

1 

C24H20 

30 (348) 

2 

10(116) 

8 

C24H21 

- 

1 

- 

1 

C24H22 

52 (603) 

1 

10(116) 

10 

C24H23 

- 

2 

- 

1 

C24H24 

33 (382) 

1 

10(116) 

12 

C24H25 

- 

1 

- 

1 

C24H26 

33(382) 

2 

10(116) 

14 

C24H27 

- 

1 

- 

1 

C24H28 

39 (452) 

1 

10(116) 

16 

C24H29 

- 

1 

- 

1 

C24H30 

34 (394) 

2 

10(116) 

18 

C24H31 

- 

1 

- 

1 

C24H32 

49 (568) 

2 

10(116) 

20 

C24H33 

- 

2 

- 

1 

C24H34 

49 (568) 

1 

10(116) 

22 

C24H35 

- 

1 

- 

1 

C24H+ 

- 

12 

- 

- 

C24HI 

- 

11 

10(116) 

1 

C24H+ 

- 

10 

10(116) 

2 

C 24 Ht 

- 

9 

10(116) 

3 

C24HI 

- 

8 

10(116) 

4 

C24H+ 

- 

7 

10(116) 

5 

C24H; 

- 

6 

10(116) 

6 

C 24 H,^ 

- 

5 

10(116) 

7 

C24HJ 

- 

4 

10(116) 

8 

C24HJ 

- 

3 

10(116) 

9 

C 24 H|o 

- 

2 

10(116) 

10 

C 24 Hr, 

- 

1 

10(116) 

11 

C24HI2 

- 

12 

10(116) 

12 

C24HI3 

28 (324) 

1 

- 

1 

C24HI4 

- 

2 

10(116) 

2 

C24HI, 

46 (533) 

1 

- 

1 

C 24 H|, 

- 

1 

10(116) 

4 

C24HI, 

64 (742) 

2 

- 

1 

C 24 H|g 

- 

1 

10(116) 

6 

C 24 Hrg 

35 (406) 

1 

- 

1 

C24H70 

- 

2 

10(116) 

8 


Article number, page 11 ofO 
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Hydrogenation 


Abstraction 


Reactant 

Barrier / meV (K) Reactive sites 

Barrier / meV (K) Reactive sites 

C 24 H+ 

30 (348) 

1 

- 

1 

C 24 H 22 

- 

1 

10(116) 

10 

C24H;J^, 

52 (603) 

2 

- 

1 

C 24 H +4 

- 

1 

10(116) 

12 

C24H+g 

33 (382) 

1 

- 

1 

C24H2-, 

- 

2 

10(116) 

14 

C 24 H 27 

33 (382) 

1 

- 

1 

C24H2g 

- 

1 

10(116) 

16 

C 24 H 29 

39 (452) 

1 

- 

1 

C 24 H 3 Q 

- 

2 

10(116) 

18 

C24Ht, 

34 (394) 

1 

- 

1 

C 24 H 32 

- 

2 

10(116) 

20 

C 24 H 33 

49 (568) 

2 

- 

1 

C 24 H 3+4 

- 

1 

10(116) 

22 

C24Ht, 

49 (568) 

1 

- 

1 

C 24 H 36 

- 

- 

10(116) 

24 

C24Ht, 

- 

- 

10(116) 

24 


Table A.2. Integrated photodissociation rates for the different hydrogenation and ionization states of the c oronene molecule. E n is the energy 
necessary for the proc ess and A5 is the ener gy d ifference between the transition state and the initial molecule iLe Page et alJl20()Ih . Values for Eo 
and A5 are taken from iLe Page et all ( 1200 Ih and iParis'et al] ( 1201411 . 


Reaction 

Eo (eV) 

AS (cal K-^) 

^ph (S 

‘G 

4*) 

Cor + hv ^ CorH_i + H 

4.73 

5.0 

9.29- 

10 - 

10 

Cor + hv ^ CorH -2 + H 2 

4.75 

5.0 

9.29- 

10 - 

10 

Cor+ + hv ^ CorH+j + H 

4.85 

5.0 

3.38- 

10 - 

10 

Cor+ + hv ^ CorH+, + H 2 

4.82 

5.0 

3.38- 

10 - 

10 

CorH_n + hv —> CorH_(n+i) + H 

3.56 

5.0 

3.49 

10 

-9 

CorH^n + hv —> CorH_(n+ 2 ) + H 2 

3.56 

5.0 

3.49 

10 

-9 

CorHl„ + hv ^ CorH!(„_^j) + H 

4.28 

5.0 

5.05 

10 

-9 

CorHl;,, + hv —> CorH^j^^ 2 ) + ^2 

4.28 

5.0 

5.05 

10 

-9 

CorHn + hv —» CorHn_i + H 

1.2 

5.0 

2.58 

10 


CorHn + hv —> CorHn -2 + H 2 

1.6 

5.0 

1.23 

10 

-9 

CorH+ + hv ^ CorH+ j + H 

2.9 

5.0 

2.52 

10 

-8 

CorH+ + hv ^ CorH:, + H 2 

3.2 

5.0 

4.91 

10 

-9 
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